OPTIMAL AMBIGUITY FUNCTIONS AND WEIL'S EXPONENTIAL 

SUM BOUND 
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Abstract. Complex-valued periodic sequences, u, constructed by Goran Bjorck, are an- 
alyzed with regard to the behavior of their discrete periodic narrow-band ambiguity func- 
tions Ap(u). The Bjorck sequences, which are defined on Z/pZ for p > 2 prime, are 
^ unimodular and have zero autocorrelation on (Z/pZ) \ {0}. These two properties give rise 

^"5 to the acronym, CAZAC, to refer to constant amplitude zero autocorrelation sequences. 

(O The bound proven is |Ap(u)| < 2/y/p + A/p outside of (0,0), and this is of optimal magni- 

^~~^ tude given the constraint that m is a CAZAC sequence. The proof requires the full power 

^^ of Weil's exponential sum bound, which, in turn, is a consequence of his proof of the Rie- 

,^^ mann hypothesis for finite fields. Such bounds are not only of mathematical interest, but 



o 



u 



X 



they have direct applications as sequences in communications and radar, as well as when 
the sequences are used as coefficients of phase-coded waveforms. 



1. Introduction 



^ 1.1. Purpose. Let Z denote the ring of integers and let C denote the field of complex 

r^ numbers. Given an integer A^, form the ring Z/A^Z of integers modulo A^. 

00 

OO Definition 1.1. Let u : Z/A^Z :— )■ C be an A^-periodic sequence. The discrete narrow band 

ambiguity function, A]\f{u) : Z/A^Z x Z/A^Z — )■ C, is defined to be 

O 1 ^-1 



AN{u)[m,n] = ^ XI "["" + k]u[k]e-'''''^''/^ 



N 

k=0 



for all (m, n) e Z/NZ x Z/NZ. 

The discrete autocorrelation of u is the function An{u)[-, 0] : Z/NZ — > C 



The ambiguity function in Definition LI stems from P. M. Woodward's definition of the 
narrow band ambiguity function defined on M x M |37] . 

Definition 1.2. An A^-periodic sequence u : Z/NZ — )■ C is constant amplitude zero auto- 
correlation {CAZAC) if it satisfies the following properties: 

(CA) Hk]\ = l ioT aWkeZ/NZ, and 

1^' 



(ZAC) C{u)[m] = —2_^u[m + k]u[k] = for all m G Z/A^Z \ {0}. 

A:=0 

Clearly, C{u)[m] = AN{u)[m,0] for each m G Z/NZ. Equation (CA) is the condition 
that u has constant amplitude 1. Equation (ZAC) is the condition that u has zero auto- 
correlation. 
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Our setting is almost exclusively limited to the case that N = p is prime. As such, Z/pZ 
is a field. 

We shall use a remarkable construction of CAZAC sequences Up of prime length p to 
prove optimal behavior of Ap{up). The construction is due to Goran Bjorck [8] (1985), 
[9j(1990). By optimal behavior, we mean that if p is an odd prime, then 

2 4 
(1) |Ap(Mp)[m,n]| < — + - for all (m,n) G (Z/pZ X Z/pZ) \ {(0,0)}, 

y/P P 



see Theorem |3.8[ By comparison, a short and elementary calculation shows that for any 
CAZAC u, 



max 



{\Ap{u)[m,n]\ : {m,n) G {Z/pZ x Z/pZ) \ {(0,0)}} > 



Vp^^' 

and therefore the bound ([I]) above is indeed of optimal order of magnitude. 



Remark 1.3. The proof of Theorem 3.8 requires Andre Weil's exponential sum bound, [35] , 
which is a consequence of his proof of the Riemann Hypothesis for curves over finite fields, 
[36], announced in the Comptes Rendus in 1940. Further, there are no more elementary 
means to prove the inequality (jTl). In fact, in estimating Ap{up), the critical term to 
estimate is a Kloosterman sum; and, if there were an easier way to bound it by C/^/p, then 
there would be an easier way to prove Weil's bound for Kloosterman sums, which is an 
essential consequence of [35] and which has withstood the test of time vis a vis evolutionary 
simplification. 

Remark 1.4. Notwithstanding the level of mathematics required to prove the inequality 
([T|, as noted in Remark 1.3, we emphasize that our coding and implementation of Bjorck's 



CAZAC sequences is truly elementary. In this regard, see [4], as well as earlier Bjorck 
experiments and constructions by one of the authors, e.g., see [7] and references therein. 



cf. Remark L5 In the parlance of waveform design. Theorem |3.8| is an ideal discrete 
"thumbtack" narrow band ambiguity function which can be used to design ideal phase- 
coded waveforms devoid of any substantial time or doppler coupling in the continuous 
narrow band ambiguity function plane. With regard to hardware implementation of these 
phase-coded waveforms (as well as others stemming from low correlation sequences), the 
power, bandwidth, and hardware requirements will introduce noise. It is understood that 
modifications must be made to the formulation of a given low correlation sequence to permit 
implementation while controlling this noise. 

1.2. Background. The study of CAZAC sequences and of other sequences related to 
optimal autocorrelation behavior has origins in several important applications, one of the 
most prominent being in the general area of waveform design associated with radar and 
communications, see, e.g. , according to year of publication [201 ED [HI E31 EH [IHl ED [23 1211 
\W\ [21 [21 [221 [121 13 [2Z]- There are hundreds of articles in this area and so this selection may 
seem arbitrary, although several of these references contain focused lists of contributions 



and specific applications. Also see Remark 1.5 



There are also purely mathematical origins for the construction of CAZAC sequences. 
One such origin is due to Norbert Wiener, e.g., see new related constructions in [HI [5]. 
Another may be said to have originated in a question by Per Enfio in 1983. This particular 
mathematical path has been documented and built upon by Bahman Saffari [21]. Enfio's 
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question is the following for a given odd prime p. Is it true that the Gaussian sequences, 
u : Z/pZ —7- C, defined by 



^rk^+sk 



u[k]=CT^-'\ A; = 0,l,...,p-1, 



where (p = e^'^*/^, r,s E Z and p does not divide r, are the only unimodular sequences of 
length p, with u[0] = 1, whose Discrete Fourier Transform (DFT) has modulus 1? This is 
equivalent to asking whether such sequences u with m[0] = 1 are the only bi- unimodular 
sequences of odd prime length. Enflo was interested in this because of a problem dealing 
with exponential sums. 

Enflo's question has a positive answer for p = 3 and p = 5. In 1984, by computer 
search, Bjorck discovered counterexamples to the Enflo question for p = 7 and p = 11, see 
[H]. Later in 1985, Bjorck saw the role of Legendre symbols in his counterexamples, and 
this led to his theorems in [9j . It also led to a host of mathematical problems, many still 
unresolved, about the number of CAZAC sequences for a given length, see, [HI El El S], as 
well as a several valuable oral and email communications by Saffari 



Remark 1.5. It is relevant to mention a striking recent application of low correlation 
sequences to radar in terms of compressed sensing [17]. In this case, the authors use Alltop 
sequences fl\ Theorem 2, cf. [32] Section 2.1.3. It then becomes natural to think in terms of 
frames generated by Bjorck sequences for extending the high-resolution radar/compressed 
sensing setting of the authors of [T7] . 



Another approach to the problem addressed in Section 1.1 is found in [13], cf. [23] ■ The 
authors obtain bounds comparable to those found herein, but their class of signals, called 
the oscillator system, is not necessarily ZAC although excellent cross-correlation criteria 
are obtained, something we have not pursued. More important, from the point of view 
of application, the characterization and construction of the oscillator system are decidedly 
representation theoretic. As such an explicit algorithm associated with the collection of 
split tori in Sp requires a Bruhat decomposition. 

The companion, ^, of this paper not only exhibits the simplicity of implementation 



stressed in Remark 1.4, but also reflects the combinatorial and geometrical complexity in 
the ambiguity function domain due to the role of the Legendre symbol in defining Bjorck 
sequences. Some of this complexity is characterized by intricate Latin and magic square 
patterns. Further, the simplicity of implementation gives rise to useful, efficient bounds off 
of small neighborhoods of (0, 0) in the ambiguity function domain for compactly supported 
waveforms on M having p lags whose coefficients are the elements of a Bjorck sequence Up. 
Also, it is not difficult to see that, as with the oscillator system, there is Fourier invariance 
of Bjorck sequences, most simply calculated in the p = 1 (mod 4) case, e.g. 



1.3. Outline. We define Bjorck sequences in Section [2j Properties of Kloosterman sums 



are proven in Section 3J_ These, in turn, are used along with Weil's results and the 
proper decomposition formula to express Bjorck sequences in the way that allows us to 
prove Theorem 3^ in Section 3^ Section |4| provides figures and data which motivated and 
guided us. 
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2. BJORCK SEQUENCES AND MULTIPLICATIVE CHARACTERS 

For each prime number p, recall that the Legendre symbol modulo p is the function 



X= (-1 :Z/pZ^ {+1,0,-1} given by 

I +1 a k = m? (mod p) for some m E Z/pZ^ , 
X[k]= l-j = lo iik = (modp), 

1—1 a k ^ m^ (mod p) for all m E Z/pZ. 

The preimage of +1 under the Legendre symbol function is the set Q of nonzero quadratic 
residues modulo p; and the preimage of —1 under the Legendre symbol function is the 
set Q^ of quadratic nonresidues modulo p. Among the many properties of the Legendre 
symbol, we shall use the fact that it is a character of the multiplicative group (Z/pZ) ^ . 
This means that Xi when restricted to (Z/pZ)^, is a group homomorphism into C^; see 
[I5j . Chapters V and VI. 



Definition 2.1. The Bjorck sequence of length p, where p is a prime and p = 1 (mod 4), 
is defined by 

u[k] = exp {i9x{k)) = exp I ^^ I — ) ) , where 9 = arccos 



for all k E Z/pZ. 

The Bjorck sequence of length p, where p is a prime and p = 3 (mod 4), is defined by 

fl — p 

a k E Q^ C (Z /pZ) ^ , where = arccos I 

\l+p 

otherwise, 




In the case p = 1 (mod 4), Definition 2.1 is equivalent to the following definition for the 
Legendre symbol sequence {0, 1, ..., —1, ..., 1} of length p. We replace the first term by 1, 
every term 1 by 

rj = exp I I arccos = \- 1- 



P-1 / Vp + 1 v^+1 

and every term —1 by the complex conjugate of r^; see |28] for a modest generalization. 
As proven by Bjorck and differently in [4J, we obtain a CAZAC, and hence bi-unimodular, 
sequence with three values, viz., 1 at fc = 0, and r] and f] ai k E (Z/pZ)^. 

In the case p = 3 (mod 4), Definition 2.1 is equivalent to the following definition for the 
Legendre symbol sequence {0,1,...,— 1} of length p. Replace the first term by 1, and 
replace every —1 by 

4 = exp t arccos = h «: 



pj I + p I + p 

As proven by Bjorck and differently in [4j, we obtain a CAZAC, and hence bi-unimodular, 
sequence with only two values, viz., 1 and ^. 
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3. The main theorem 
3.1. The Legendre symbol and Kloosterman sums. 
Definition 3.1. Let p be a prime. For any integers a, b, ttie quantity 

K[a,b;p]= >^ exp (2TTi{ax + bx^^)/p), 

where x~^ denotes the multiplicative inverse of x in the field Z/pZ, is a Kloosterman sum. 

Kloosterman sums are always real-valued, as the following Lemma states. 
Lemma 3.2. Let p be a prim,e. Then K[a, b\p\ G M for all integers a, 6 G Z. 
Proof. By the substitution y = —x, we have 

K[a,b;p] = ^ ^~2ni{ax+bx-^)/p ^ J2 e^^'^^y+^y^'^IP = K[aMp\- □ 

The following classical description of certain Kloosterman sums was first observed by 
Hans Salie in equation (52) of [30], using a formula of Ernst Jacobsthal from a footnote 
on page 239 of yjj. Jacobsthal's footnote refers the reader to his 1906 Ph.D. thesis, but 
fortunately the proof of his formula is not difficult to derive. 



Lemma 3.3. Fix an odd prime p and an integer a not divisible by p. Let x = - denote 

the Legendre symbol modulo p. 

a. (Jacobsthal, 1907) Let F : Z/pZ -^ C be any function. Then 

p-i p-i 

xeZ/pZ'x x=0 x=0 

p-1 

b. (Salie, 1932) K[l,a;p] = ^ x[a;^ - 4a]e2"*^/P. 

x=0 



The formulas of Lemma 3.3 are known, but we include their proofs because of the role 
they play in our approach. 

Proof, (a). Let g : (Z/pZ)^ — )• Z/pZ be the function g[x] = x + ax~^. For each t G Z/pZ, 
set 

N[t] = card{x G {Z/pZ^ : g[x] = t}. 
The desired sum can now be written as 

J2 F[x + ax-']= J2 F[g[x]]= Yl ^M^W- 
xez/pzx xez/pzx tez/pZ 

Thus, it suffices to show that N[t] = 1 + xl^"^ ~ ^o]. 

Note that g[x] = g[ax^^]. Conversely, for any x, y G (Z/pZ)^ with g[x] = g[y], we must 
have either y = x or y = ax~^, since = g[x] — g[y] = {x — y){l — ax~^y~^). Thus, 
N[t\ < 2 for all t G Z/pZ, and N[t] = 1 if and only if t = g[x] for a point x G (Z/pZ)^ such 
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that X = ax~^. This latter condition occurs if and only ii x^ = a in Z/pZ; in that case, 
t = g[x] = g[ax'^] = 2x, or equivalently, t^ = 4a. Thus, if we set 

S = {g[x] : X e (Z/pZ)^a;2 ^ a}, 

then 

(2 iiteS, 

N[t] = ll if t2 = 4a, 

I otherwise. 
Note, on the other hand, that 

{2 if t^ — 4a is a square in (Z/pZ)^ , 
1 if t^ _ 4a = in Z/pZ, 
otherwise. 
Thus, it suffices to show that 

(2) S = {te Z/pZ : t^ - 4a is a square in {Z/pZ^}. 

Given t G 5, pick x G (Z/pZ)^ such that t = g[x]. Then 

t^ — Aa = ix + ax^'^Y — 4a = x^ — 2a + c?x~'^ = (x — ax~^)^. 



X 



In addition, since t^ 7^ 4a for all t G S*, it follows that t^ — 4a is a square in (Z/pZ) 
proving the forward inclusion. 

Conversely, given t G (Z/pZ) for which there is some z G (Z/pZ)^ with z"^ = t"^ — 4a, 
set X = {t + z)/2 G (Z/pZ)^. Then x(x — z) = (t^ — z'^)/A = a, and therefore g{x) = 
X + ax~^ = 2x — z = t. It follows that t E S, proving equation pi) and hence part (a). 

Part (b) is immediate by setting F[x] = e^'^^^lp and noting that X]x=o e^'^*^/^ = 0. D 



Theorem 3.4. Fix an odd prime p. Let x = { ~ ] be the Legendre symbol modulo p. Then 

e-™"/^Ap(x)[m,n] G M, and \AM[m,n]\ < -^, 

for all m,n E Z/pZ \ {0}. 

Proof. Fix m,n E Z/pZ \ {0}. Noting that x is multiplicative and real-valued, we have 



AM[m,n] = - J2 x[A: + m]x[fc]e-2-^^"/^ = - ^ x[Hk + m)]e- 
^ fcez/pZ ^ fcez/pZ 



2-iTikn/p 



Let a = (r7in)^/16, 6 = m/2, and c = — 1/n, where we are doing the arithmetic in Z/pZ. 
Substituting k = ex — b, we have 

y4p(x)[m,?2] = - y^ x[(cx — b){cx + b)~\ exp(— 27rm(cx — b)/p) 
^ xez/pZ 

„27rifen/p p27rjbn/p 

(3) = Yl ^[^'^' - ^']^'^'^^' = ^ ^[i,«;p], 



where the final equality is valid because 6^ = 4ac^, and hence 

X[c^x^ - b^] = x[c\x^ - 4a)] = x[c^]x[x^ - M = x[x^ - 4a]. 
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Since [e^^ibn/py ^ ^2nimn/p ^ (^^nimn/p^2^ ^^ j^^^g ^2nibn/p ^ _^^nimn/p ^ ^^^ therefore by 
equatfon ^ and Lemma 3.2 



g-TT*^ 



'^'""/Mp(x)[m,n] = ±-K[l,a;p] G M. 



Finally, because a G Z/pZ \ {0}, we have |i^[l,a;p]| < 2y/p, by Weil's bound for Kloost- 
erman sums in [35j. Thus, equation ^ gives |Ap(x)[m, n]! < 2/^, as desired. D 



Remark 3.5. In [35], Weil proves his bound for |K[a,6;p]| by first using Lemma 3.3 



to 

rewrite K[a,b;p] as '^xl^'^ ~ 4a]e^'^*'^'/^ and then bounding the new sum. Philosophically, 
then, it would be more direct not to convert the sum to the form ^exp(27ri(x + ax~^)/p). 
Nevertheless, we have applied the transformation in Lemma 3^ because the latter form of 
Kloosterman sums is better known than are the details of Weil's proof. 

3.2. Main bound. We shall need the following technical lemma, which gives bounds for 
the ambiguity function of any sequence that is a function of the Legendre symbol. 

Lemma 3.6. Fix an odd prime p and complex numbers r,s,t G C. Let x '■ Z/pZ — )■ C 6e 
the Legendre symbol modulo p, and let U : Z/pZ — )■ C 6e the function 

{r ifx{k) = l, 
s z/x(A:) = -l, 
t ifk = 0. 

Set R={r + s)/2, S = {r - s)/2, T = t - R, and Cp = e^^'/P. Then 

Ap{U)[m,n] = \S\'^Ap{x)[m,n] + -[Ei[m,n] + E2[m,n]) 

for all m,n E TLjpTL \ {0}, where -E'i[m, n] = RT + RTCJ^'^ , and 

E\mn] = [^^^ + ^TCr)xH + {RS + RSC'^^)x[n]^ tf p = 1 (mod 4), 
2[m,nj y^sf-STC;''')x[m]-{RS + RSC'^^)tx[n]^ zf p = 3 (mod 4). 

Proof. For any two functions F, G : Z/pZ — t- C, write 

BJF, G) [m, n] = - V F[k + mlGMe-^"^'^"/^. 
^ kei/pi 

Define functions rj, 5 : "L/pl, — )■ C by 

Thus, U = Rrj + Sx + T6, and hence 

Bp{U,U) = \R\'Bp{ri,r]) + \S\'Bpix,x) + \T\'Bp{S,S) 
+ RfBpiri, 5) + RTBpiS, r/) + SfBp{x, 5) + STBp{5, x) + RSBp{r], x) + RSBp{x, v)- 
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To compute Bp{U, U), we shall compute each of these nine terms separately. Since m ^ 0, 
we have Bp{6, 6) = 0. In addition, Bp{ri, rj) = 0, since J2k(^z/pz e"^'^*'^"/^ = and n 7^ 0. We 
also have Bp{x,x) = ^p(x) by definition. Meanwhile, it is immediate that 

pBp{r], 6) [m, n] = 1, pBp{6, r]) [m, n] = Cp™", 

pBpix,S)[m,n] = x[m], and pBp{6,x)['m',n\ = C™"x[-"^]- 

Next, p-Bp(?7, x)[^,''^] = Ti—n'jp], where r[a;p] is the Gauss sum 

T[a-p]= J2 xWe'™'^/^ 
kez/pZ 

However, Gauss proved that r[a;p] = sx['^]y/P^ where e = 1 ii p = 1 (mod 4), and e = i 
if p = 3 (mod 4); see, for example. Proposition 6.3.1 and Theorem 6.4.1 of [18]. Hence, 

pBp{v,x)[m,n] =ex[-n]y/p. Similarly, 

pBp{x,v)[^,n]= J2 X[k + m]e-'-'''^/P = J^ ^[j]e'^-^U-m)n/p 
kez/pZ j&^/p^ 

= Cr^[-n;p]=eCrx[-n]VP- 
Combining the nine computations above, and noting that 

r ,1 r Ti r,iJxM ifp^l (mod 4), 
I— X[A;J np = 6 (mod 4), 

we have Bp(U, U) = \S\'^Ap{x) + (-E'l + E2)/p, where Ei and E2 are the quantities in the 
statement of Lemma 13.61 D 

The following elementary bound will be needed to prove the p = 3 (mod 4) case of 
Theorem 13.81 



Lemma 3.7. Let X,Y eM., and let z E C with \z\ = 1. Then 

\zX + (1 - z^)Y\ < VX2 + 41^2. 
Proof. Noting that 2;^ = 1, we have 



zX + {1 - z^)Y\ = J{zX + {l-z^)Y){zX + {l-z^)Y) 



X2 + (2(1 - Z^) + Z{1 - Z^))XY + (1 - Z2)(l _ ^2)y2 



= ^/X^ + \l-zYY^< VX^ + AY^, 

since z{l — z"^) + ^(1 — z"^) = z — z + z — z = and |1 — z^| < 2. D 

We are now ready to state and prove our main result. 

Theorem 3.8. Let p be an odd prime, and let Up he the Bjorck function for p. Then the 
ambiguity function, Ap{up), defined on TLjpTL x Z/pZ as 

Ap{up)\m, n] = - ^ Up\k + m]^e-2"''="/^ 

^ k&LlpZ 
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satisfies the estimate 



\A{up)[m,n]\<^+^P^ 



.p 



if p = 1 (mod 4), 
^3/2 ^fP = ^ (mod 4), 



for all {m,n) E {Z/pZ x Z/pZ) \ {(0,0)}. 

Proof. Fix {m,n) G iZjpZ x Z/pZ) \ {(0, 0)}. If m, = 0, then n 7^ 0, and we have 



P 



k& 



P 



smce \u 



pl/ljI = 1 for all k G Z/pZ. On the other hand, if n = 0, then ?ti 7^ 0, and we have 
Ap[up)[ni, 0] = - y, ^p[^ + '^I'Upi^] = 0, 



P 



k&Z/pl 



because Up has zero autocorrelation. Thus, by the fact that Up is a CAZAC, we may assume 
for the remainder of the proof that m,n ^ 0. 

If p = 1 (mod 4), then in the notation of Lemma 3.6, we have r = (1 + y/p)~^{l + 

Vp 



i^2^ + p), s = (1 + y/p)-\l - i^2^ + p), and t = 1. Thus, 

1 



R 



r + s 



i + Vp' 



s 



r — s 



p + p 



1 + v^ 



and T = t-R 



1 + v/p' 



The quantities Ei and E2 in Lemma 3^ are therefore 



El [m, n] 



P I 



(1 + ^f 



and 

E2\m,n\ 



;i + vp)^ 
^fp 



1 - Cr)v^V2v^ + P- '^[^] + C - ^)M'i^^p-^x\n, 



^(i-Cr)UH-xM)V2v^ + p 



(1 + ^Y 

Noting that |1 + C™"|, |1 — C™"l; ^'^'^ Ixl"^] ~ xM I are each less than or equal to 2 and that 
^/^/p~+P < a/1 + '^y/p + P = 1 + v^P) we obtai 



am 



i£a™.„i +^.i™,„ii < i^ + ^ < (tI^ +4. 



Hence, by Lemma 3^ and Theorem |3.4[ we have 

2 2 4 

' "^ "^^ ^^-vp vpi'^ + vpy p vpi'^ + vpy 



2 4 

{2,/^ + p+l) + - 



^(1 + ^y 



:i + vp)' + 



+ 



4 



p \/p p 
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Similarly, if p = 3 (mod 4), then r = 1, s = {1 + p) ^(1 — p + 2iy/p), and t = 1, and 
therefore 



R 



r + s 



1 



S 



r — s 



-^Vp 



2 l-iy/p 2 1-iy/p 

Thus, the quantities Ei and E2 in Lemma [3. 6| are 



and T = t — R 



-iyp 
1 -ly/p' 



Ei[m, n] 



p + l 



and 

E2[m,n 



1 



p+l 
Since jS'p = p/{p + 1), we have 

\S\^Ap{x)[m,n] + -E2[m,n] 



{p - pcr)xH - i^vp - cr'VpMn]Vp 



p(i - c; 



p I 



p+l 



X\m\ + xM 



p 



p + l 



pA,{x){mA + (1 - Cr)(^N + XN) 



Setting z = e^™"/^, X = e" 



TTimn/p, 



'pylp(x)[m,n], and F = x[m] + xN? so that X G M with 



-^1 < 2y^ by Theorem 3.4, F G M with \Y\ < 2, and \z\ = 1, Lemma 3.7 tells us that 



\S\''Ap{x)[m,n] + -E2[m,n] 



P 



^ yx^ + AY^ ^ yAp + 16 _ 2yj^+4: 



p+l 



p+l 



p+l 



Hence, by Lemma 3.6 and the fact that |1 — CI""| < 2, we obtain 



\Ap{up)['m,n]\ < \S\'^Ap{x)[m,n] + -E2[m,n] 



P 



-Eilrrijn] 



^2Vp + A^ 



p + l yp{p + 1) ^(p + 



— (Vp2 + 4p+l) 



^ 2(p + 3) _ 2 ^ 



2 4 

< ^ + 



v^(p + i) Vp Vp(p + i) Vp f 



,3/2' 



D 



Remark 3.9. The bounds in Theorem |3.8| may be improved very slightly but at the 

great expense of simplicity. For example, if p = 1 (mod 4), then the bounds |1 — C™"! < 

2 and |1 + CJ^"! ^ 2 could be improved, as obviously these quantities cannot both be 

simultaneously close to 2. However, the resulting bound is far more complicated to write, 

and the savings is only about 2p~^/^, as illustrated by considering (I^"' very close to —1. 

Similarly, removing the simplification 4^^/(1 + ^/p) < 4 would also only save us about 
4p-3/2_ 



4. Figures and table 



Natural algebraic and analytic calculations convinced us that the proof of Theorem 3.8 
depended on substantial number theoretic results. In parallel. Figure [T] supported the truth 
of Theorem 3.8 before we proved it. The x-axis lists the primes between 1 and 1000. The 
y axis lists the values, 

(4) max \AJur,)\m,n\\. 

(m,n)^(0,0) ^' -^'^ 
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p prime and 1<p<1000 
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Figure 1. p and raax{\ Ap{up)[m,n]\ : {m,n) ^ (0,0)} 
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Figure 2. Ap{up) for p = 13 



Figure[T]also displays the curves y = 2/y/p and y = 2/y/p+A/p for comparison. Figure[2| for 
the case p = 13, illustrates the symmetries inherent in the function Ap{up) on Z/pZ x Z/pZ. 
These are fully explained for all p in [4j; and they led to the realization of the complexity 
involved in proving Theorem |3.8[ as well as to a host of geometrical and combinatorial 
phenomena and problems. Figure |3] illustrates Theorem 3^ for the case p = 503. 

Table [I] indicates some of the finer behavior of the quantity ^, over three different 
ranges of primes. This data suggested to us that 2/y/p was very nearly the upper bound 
for |Ap(Mp)[m,r2]|, {m,n) ^ (0,0), and it helped lead us to the proof that 2/y/p + A/p 
is an upper bound. In addition, although a number of primes p = 1 (mod 4) require 
a bound larger than 2/y/p, we noted that only very few primes p = 3 (mod 4) allowed 
|Ap(tip)[m, njl > 2/y/p for {m,n) ^ (0,0). For example, p = 139 is the only such prime 
in Table [ij Our broader calculations for other primes showed that the only such primes 
between 1000 and 5000 are 1259, 2111, and 3511; the only ones between 10000 and 24360 
are 13879, 16091 and 23719; and there are none between 100000 and 105000. Moreover, for 
all seven of those primes, the maximum value of |y4p(Mp)[m,n]| — 2/y/p for {m,n) ^ (0,0) 



12 



JOHN J. BENEDETTO, ROBERT L. BENEDETTO, AND JOSEPH T. WOODWORTH 



|/\,(U503)[m.nll 
Max away from [0,0] : 0.088064 




200 
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200 
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Figure 3. Ap{up) 



for p = 503 



is still far smaller than A/p, a fact which ultimately led us to the sharper bound for p = 3 



(mod 4) in Theorem 3.8 
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